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Abstract 

The exclusive rare radiative B meson decays to orbitally excited axial- vector 
mesons K^(1270), i^i(1400) and to the tensor meson 7^1(1430) are investi- 
gated in the framework of the relativistic quark model based on the quasipo- 
tential approach in quantum field theory. These decays are considered without 
employing the heavy quark expansion for the s quark. Instead the s quark 
is treated to be light and the expansion in inverse powers of the large recoil 
momentum of the final K** meson is used to simplify calculations. It is found 
that the ratio of the branching fractions of rare radiative B decays to axial 
vector ivrj*(1270) and -ftri(1400) mesons is significantly influenced by relativis- 
tic effects. The obtained results for B decays to the tensor meson ir|(1430) 
agree with recent experimental data from CLEO. 

PACS number(s): 13.20.He, 12.39.Ki 

I. INTRODUCTION 



Rare radiative decays of B mesons represent an important test of the standard model 
of electroweak interactions. These transitions are induced by flavour changing neutral 
currents and thus they are sensitive probes of new physics beyond the standard model. 
Such decays are governed by one-loop (penguin) diagrams with the main contribution 
from a virtual top quark and a W boson. Therefore, they provide valuable information 
about the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vts and Vtb- The statis- 
tics of rare radiative B decays considerably increased since the first observation of the 
B K*'~f decay in 1993 by CLEO |jl[. This allowed a significantly more precise deter- 
mination of exclusive and inclusive branching fractions 0. Recently the first observation 
of the rare B decays to the orbitally excited strange mesons has been reported by CLEO 
0. The branching fraction for the decay to the tensor i^'2(1430) meson has been measured 
BR{B ir*(1430)7) = {l.QQtoil ± 0.13) x 10-^ as well as the ratio of exclusive branch- 
ing fractions r = BR{B K;{U30)-f)/ BR{B K*{892)-f) = 0.39tg:i^. The data for the 
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other decay channels will be available soon. This significant experimental progress provides 
a challenge to the theory. Many theoretical approaches have been employed to predict the 
exclusive B — >■ i^* (892)7 decay rate (for a review see 0] and references therein). Consid- 
erably less attention has been payed to rare radiative B decays to excited strange mesons 
Most of these theoretical approaches [^,0 rely on the heavy quark limit both for the 
initial b and final s quarks and the nonrelativistic quark model. However, the two predic- 
tions [^,0 for the ratio r differ by an order of magnitude, due to a different treatment of the 
long distance effects and, as a result, a different determination of corresponding Isgur-Wise 
functions. Only the prediction of Ref. [[^ is consistent with the available data. Nevertheless, 
it is necessary to point out that the s quark in the final K* meson is not heavy enough, com- 
pared to the A parameter, which determines the scale of I/ttlq corrections in heavy quark 
effective theory p|. Thus the l/m^ expansion is not appropriate. Notwithstanding, the ideas 
of heavy quark expansion can be applied to the exclusive B K*{K**)'y decays. From the 
kinematical analysis it follows that the final K*{K**) meson bears a large relativistic recoil 
momentum |A| of order of mb/2 and an energy of the same order. So it is possible to ex- 
pand the matrix element of the effective Hamiltonian both in inverse powers of the b quark 
mass for the initial state and in inverse powers of the recoil momentum |A| for the final 
state P,pi][. Such an expansion has been realized by us for the B —>■ K*{892)'j decay in the 



framework of the relativistic quark model 0. In Refs. |10] it was shown that in the leading 
order of this expansion a specific symmetry emerges which imposes several relations between 
the form factors of semileptonic and rare radiative B decays. It is important to note that 
rare radiative decays of B mesons require a completely relativistic treatment, because the 
recoil momentum of the final meson is large compared to the s quark mass. The calculated 
branching fraction for this decay was found in good agreement with experimental data. 
In Ref. [|TT|] we considered the exclusive rare B decay to the orbitally excited tensor meson 
i^^^(1430). Here we extend this analysis to the axial- vector mesons i^*(1270) and i^'i(1400). 

Our relativistic quark model is based on the quasipotential approach in quantum field 
theory with a specific choice of the quark-antiquark interaction potential. It provides a 
consistent scheme for the calculation of all relativistic corrections at a given f^/c^ order and 
allows for the heavy quark l/mq expansion. In preceding papers we applied this model to 
the calculation of the mass spectra of orbitally and radially excited states of heavy-light 
mesons , as well as to the description of weak decays of B mesons to ground state heavy 



and light mesons |T3|,|Tj]. The heavy quark expansion for the heavy-to-heavy semileptonic 
transitions |TB|,|TB[ was found to be in agreement with model-independent predictions of the 
heavy quark effective theory (HQET). 

The paper is organized as follows. In Sec. II we define the form factors, which govern 
the exclusive rare radiative B decays to orbitally excited K mesons. The relativistic quark 
model is described in Sec. HI, and in Sec. IV it is applied for the calculation of the rare 
radiative decay form factors. Our numerical results for the form factors and decay rates as 
well as comparison of these results with other theoretical predictions and experimental data 
are presented in Sec. V. There we also discuss the relations between the form factors of rare 
radiative and semileptonic B decays in the formal heavy quark limit. Sec. VI contains our 
conclusions. 
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II. RARE RADIATIVE B DECAYS 



In the standard model, B decays are described by the effective Hamiltonian, after inte- 
grating out the top quark and W boson and using the Wilson expansion |T^. For the case 
of 6 ^ s transition: 

AC ^ 

Heffib ^s) = -^V*Vt,j:C,if^)OM)^ (1) 

where Vij are the corresponding CKM matrix elements, {0^} are a complete set of renor- 
malized dimension six operators involving light fields, which govern b —>■ s transitions. They 
include six four-quark operators Oj {j = 1,...,6), which determine the non-leptonic B 
decay rates, the electromagnetic dipole operator 

O7 = j^sa^'imtPR + m,PL)bF^u, Pr,l = (1 ± 75)/2, (2) 

and the chromomagnetic dipole operator Og- O7 and Og are responsible for the rare B 
decays b ^ S'j and b sg, respectively Ijl^. The Wilson coefficients Cj(/i) are evaluated 
perturbatively at the W scale and then they are evolved down to the renormalization scale 
fi ~ rrih by the renormalization group equations. There are also power-suppressed terms 

The dominant contribution to B —>■ K**'j decay rates come from the electromagnetic 
dipole operator O7. The matrix elements of this operator between the initial B meson 
state and the final state of the orbitally excited K** meson have the following covariant 
decomposition 



{K,ip\e)\MKa,,b\Bip)) = Uik') ({e* ■ k){p + p'), - tlip" - p" 

+to{e){e* ■ k){{p'-p")k, - {p + p'),e)/{MBMK,), 
{Ki{p',e)\sik,a^,-f,b\B{p)) = iU{e)e^,x.e*^k^{p + p'f , (3) 

{Kl{p\e)\sik,a,,b\B{p)) = s+ik') ((e* ■ k){p + p'), - e;{p' - p")) 

+s.{e) {{e*-k)k,-e;e) 

+soik')ie* ■ k)iip'-p")k, - {p + p'),e)/iMBMK*), 
{Kl{p',e)\mk,a^,j,b\B{p)) = is+{e)e^,xae*''k\p + p'y , (4) 

{K;{p\e)\7stk,a,,b\Bip)) = ig^{k^)e^,x^e*^^ ^k\p + p'Y , 

Mb 

{K;{p', e)\Mk.a,.j,b\B{p)) = g^i^) (^.^(p + p'), - e^^^y - p")^ 



+9o{mp' - p")k, - iP + P'),k')e*,^J^Jf-^^ (5) 

B Ki 



where e^{e^y) is a polarization vector (tensor) of the final axial- vector (tensor) meson and 
k = p — p' is the four momentum of the emitted photon. The exclusive decay rates for the 
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emission of a real photon (fc^ = 0) are determined by form factors t+(0), ^4.(0) and g+{0). 
They are given by 

= ^G-,n^|V^.V'.,nC.(,„.)r«i(0,^ (1 - (1 + . (8) 

where C^{mb) is the Wilson coefficient in front of the operator O7. It is convenient to 
consider the ratio of exclusive to inclusive branching fractions, for which we have 



^ ^ BR{B iri(1400)7) _ (1 - MIJMI) (l + Mj^Mj^ 

, 3 



^ BR{B ir-(1430)7) 1 , Ml (l - I Ml) (l + Ml, I Ml ^ 
""'^ BR{B-^Xsi) 8^^^ 'mi, {l-ml/mlf{l + ml/ml) ' ^ ' 



The recent experimental value for the inclusive decay branching fraction PU 

BR{B Xsi) = (3.15 ± 0.35 ± 0.32 ± 0.26) x 10"^ 
is in a good agreement with theoretical calculations (see e.g. P). 

III. RELATIVISTIC QUARK MODEL 

Now we use the relativistic quark model for the calculation of the form factors t+(0), 
s+(0) and g+{0). In our model a meson is described by the wave function of the bound 



quark-antiquark state, which satisfies the quasipotential equation |21| of the Schrodinger 
type in the center-of-mass frame: 



*m(p) = / 7^V{p, q; M)^M{q), (12) 



V 2/iH 2fiRj J {2tt) 

where the relativistic reduced mass is 

_ - {ml - ml? 
^« 4AP ' ^^^> 

and }P'{M) denotes the on-mass-shell relative momentum squared 
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The kernel V{p, q; M) in Eq. (0) is the quasipotential operator of the quark- ant iquark 
interaction. It is constructed with the help of the off-mass- shell scattering amplitude, pro- 
jected onto the positive energy states. An important role in this construction is played by 
the Lorentz-structure of the confining quark- ant iquark interaction in the meson. In con- 
structing the quasipotential of the quark-antiquark interaction we have assumed that the 
effective interaction is the sum of the usual one-gluon exchange term and the mixture of 
vector and scalar linear confining potentials. The quasipotential is then defined by 



with 



^(P, q; M) = Ug{p)uQ{-p)V{p, q; M)ug{q)uQ{-q), 



V(p, q; M) = -a,D,Mli;ih + V^U^)T^,^q,, + (k), 



(15) 



where is the QCD coupling constant, D^^ is the gluon propagator in the Coulomb gauge 
and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 



e{p) + m 



( 



1 



X 



(16) 



with e(p) = -^/p^ + m^. The effective long-range vector vertex is given by 



r^(k) = 7^ + 'i^<y^uk\ 



(17) 



where k is the Pauli interaction constant characterizing the nonperturbative anomalous chro- 
momagnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 



^ { 

conf \ 



;i-5)(Ar + i?), 



conf 



(18) 



reproducing 



conf I 



Konf(0 + K 



^ { 
conf \ 



Ar + B, 



(19) 



where e is the mixing coefficient. 

The quasipotential for the heavy quarkonia. 



expanded in v^jc^^ can be found 



m 



Refs. P^IP^ and for heavy-light mesons in |T^. All the parameters of our model, such as 
quark masses, parameters of the linear confining potential, mixing coefficient e and anoma- 
lous chromomagnetic quark moment k, were fixed from the analysis of heavy quarkonia 

4.88 GeV, 



masses 



23] and radiative decays |25 



The quark masses m\, = 4.88 GeV, = 1.55 GeV, 
rus = 0.50 GeV, mu,d = 0.33 GeV and the parameters of the linear potential A = 0.18 GeV^ 
and B = —0.30 GeV have the usual quark model values. In Ref. [|T5| we have considered 
the expansion of the matrix elements of weak heavy quark currents between pseudoscalar 
and vector meson ground states up to the second order in inverse powers of the heavy quark 
masses. It has been found that the general structure of the leading, first, and second order 
l/rriQ corrections in our relativistic model is in accord with the predictions of HQET. The 
heavy quark symmetry and QCD impose rigid constraints on the parameters of the long- 



range potential in our model. The analysis of the first order corrections fixes the value of 
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the Pauli interaction constant k = — 1. The same value of k was found previously from the 
fine splitting of heavy quarkonia ^Pj- states [^]. The value of the parameter characterizing 
the mixing of vector and scalar confining potentials, e = — 1, was found from the analysis 
of the second order corrections [jl5l. This value is very close to the one determined from 



considering radiative decays of heavy quarkonia [ 25 



IV. RARE RADIATIVE B K**^ DECAY FORM FACTORS 

In the quasipotential approach, the matrix element of the weak current = s^k^a^u{l + 
7^)6 between the states of a -B meson and an orbitally excited K** meson has the form |2T 



{K**\J,mB) = I ^^^'/...(p)r,(p,q)vl/^(q), (20) 

where r^(p, q) is the two-particle vertex function and '^b,k** are the meson wave functions 
projected onto the positive energy states of quarks and boosted to the moving reference 
frame. The contributions to T come from Figs. 1 and 2. The contribution T^"^^ is the 
consequence of the projection onto the positive-energy states. Note that the form of the 
relativistic corrections resulting from the vertex function r^^-* explicitly depends on the 
Lorentz structure of the gg-interaction. The vertex functions look like 



and 



r«(p, q) = n.(pi)^a^.r(l + j')uM{2nf6{p2 - ^2), (21) 



r?np,q) = ^.(pi)^.(P2)^{^^i,.fc.(i + 7i')^rTi^4^7?v(p2 - q2) 

+^(P2 - q2) ^ + (22) 

where ki = pi - A; k'^ = qi + A; A = p^** - Pb] 

^(_) ^ e{p) - (m7° + 7°(7p)) 
2e(p) 

The wave functions of P-wave K** mesons at rest can be parametrized either through the 
wave functions of the states ^Pi, ^Po,i,2 used in quark models for quarkonia (LS-coupling 
scheme) or through the wave functions K{l/2) and K{3/2) used in HQET (js-coupling 
scheme). The structure of the wave functions for the states with 0"*" and 2+ quantum numbers 
is the same in both parametrizations, while two real states with l"*" quantum numbers are 
different mixtures of states in these parametrization. Experiment shows that iri(1400) and 
iir^(1270) mesons are nearly equal mixtures of ^Pi and ^Pi quark model states As a 
result the HQET parametrization turns out to be more appropriate since the real ii'i(1400) 
and i^^(1270) mesons almost coincide with the corresponding states in js-coupling scheme. 
The wave functions at rest in HQET parametrization are given by 

^K-{P) = KUp) = 3^rVi.0)(p), (23) 
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where J and M are the meson total angular momentum and its projection, while j is the 
M, d-quark angular momentum; ipK{j) (p) is the radial part of the wave function. The spin- 
angular momentum part 3^/^ has the following form 

xY,''-''^-''^XQ{crQ)xM. (24) 

Here {jirrii, j2m2\J M) are the Clebsch-Gordan coefficients, yj*" are the spherical harmon- 
ics, and (where a = ±1/2) are spin wave functions : 

x{im = (^X x(-i/2) = f?V 



Then 



K.nm = ^K%/2) cos + *^*fi/2) sin , 



^^'*(i270) = COS - ^],%/2) Sin 0, (25) 

where is a small mixing angle. We have calculated the wave functions of orbitally excited 
K** mesons in our model by the numerical solution of Eq. (|T^ with the quasipotential (|T3p 
expanded in inverse powers of quark energies. Such an expansion is more adequate for K 
mesons than the usual nonrelativistic expansion since the velocities of light u, d, s quarks are 
highly relativistic. The calculated spin-averaged P-wave K meson masses as well as spin- 
orbit splittings are consistent with experimental values. The obtained value of the mixing 
angle in ( P^D also agrees with experiment and is approximately equal to ~ 4°. In the 
following we will use the functions (^) for decay form factor calculations assuming that the 
physical form factors for B — > K[*'^'y decays are related to the calculated ones (denoted by 
a tilde) by 

t+ = t+ cos + s+ sin 0, 

s+ = s+ cos — t+ sin 0. (26) 

It is important to note that the wave functions entering the weak current matrix element 
(pO|) cannot be both in the rest frame. In the B meson rest frame, the K** meson is moving 
with the recoil momentum A. The wave function of the moving K** meson "^k** a is 
connected with the K** wave function in the rest frame o = "^k** by the transformation 



*x"a(p) = Dl/\RYjDl/\RYj^K,,o{p), (27) 

where is the Wigner rotation, is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D^^'^{R) in spinor representation is given by 

{l ^^Dl/J{RYJ = S-\ps,,)S{A)S{p), (28) 
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where 



i{p) + m 
2m 



1 + 



ctp 



e{p) + m J 

is the usual Lorentz transformation matrix of the four-spinor. 

We substitute the vertex functions r^^-* and T^'^^ given by Eqs. (^) and ( p2D in the decay 
matrix element (^OD and take into account the wave function transformation (^^. The 
resulting structure of this matrix element is rather complicated, because it is necessary to 
integrate both over d^p and d^q. The 6 function in expression (^) permits to perform one 
of these integrations and thus this contribution can be easily calculated. The calculation 
of the vertex function F^^^ contribution is more difficult. Here, instead of a 5 function, we 
have a complicated structure, containing the qq interaction operator V. However, we can 
expand this contribution in the inverse powers of the heavy b quark mass and large recoil 
momentum |A| ~ mb/2 of the final K** meson. Such an expansion is carried out up to 
the second order. ^ Then we use the quasipotential equation in order to perform one of 
the integrations in the current matrix element. As a result we get for the form factors the 
following expressions with k = — 1 



t+(0) = t«(0) + (1 - e)tf''(0) + 5tt^^(0), 

1 l EKi3/2) |A| r d^p - / 2eg 

3v^V Mb EKi3/2)+MKi3/2)J {27Cy ^^'^^^^^ EKi3/2) + MKi3/2) 



(29) 



ts{p + A) + m 
Mb-E 



-.- + mb / ^(P■A) 

W{3/2) 1 , [ P 



pA^ 



s(p + A) 



eg(p) + m 



^ _^ Mb — Ek(3/2) 



P 

g es{p + A)+ms 



es{p + A) + rris [e^ (p + A) + m^] [e;, (p) + nib] 

Mb + Mk(3/2) P \ , . 

Mb - Mk(3/2) esip + A) + /^^^^P^' 

1 I Ek{3/2) |A| r _fp_- / 

3v /2V Mb E k^,/,) + M^(3/2) J (2vr)3^^'(^/^HP + 

3 (-Ex (3/2) + ^iC(3/2) 



(30) 



2e. 



ea(A) + J 
^ 2e.(A) 1 
eg(A) - mg 
■2e,(A)[e,(A) + m, 
2e„ 



(p. A) 



Ek{3/2) + ^^X(3/2) 



pA^ 



(p) + 



es(A) + 
2e„ 



-es(p + 



M 



i^(3/2) 



'^(3/2) + ^^-(3/2) 



— e 



<(p + 



E 



'^(3/2) + ^X(3/2) 



^This means that in expressions for i, s, (0) and s, ' (0) we neglect terms proportional to 
the third order product of small binding energies and ratios p^/eg(A), p-^/e^(A) as well as higher 
order terms. 
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P 



eq{p) + ruq 
-es(p + 



MB-eb{p)-eq{p) es{A)-ms 



eb{A) + nib 2e,( A) [e, (A) + m. 



K{3/2) 



2e„ 



A 



a) -e,(p + 



2e„ 



Ek{Z/2) + MK{-i/2) 



^))\]Mp), (31) 



r(2)5.„x ^ J_ / £^K(3/2) |A| r (Pp - / 2eg ^ 

+ 3^2 V £^;f(3/2)+Mx(3/2)i (27r)3'^^('/')vP £;^(3/2)+M^(3/2) ^ 



X 



eg (A) +mg 
^ 26, (A) 

e,(A) - nis 
2e,(A)[e,(A)+m, 



-3(-E'i^(3/2) + ^j^(3/2)) ^ n ' + 



P 



pA2 eg(p) + mq 



Mb - eb{p) - tq{p) - 



2e„ 



-a) -e,(p + 



Mb — Ek{3/2) 
es{A) +ms 
2e„ 



M, 



K{3/2) 



Ek{3/2) + Mk{3/2) 

5+(0) = g;) (0) + (l - ^(0) + £5f ^(0), 



Ek{3/2) + Mk{3/2) 



1 E 



K{l/2) 



i'K(l/2) (P + 



2e„ 



3V Mb EKii/2) + MKii/2) J (27r)3 

3(-Ex(l/2) + Mk{1/2)) 



X 



X 



es(p + A) +ms 



\\ 2e,(p + A) ^ 2efc(p) 
Mb — Ek(i/2) 



€b{p) + mt, 



Ek{1/2) + Mk{\/2) 

P-A) 



>^b(p), (32) 
(33) 



1 + 



es{p + A)+ms 



2p 



2p 



X 



X 



^ ^ Mb - Ek(i/2) ^ 



eq{p)+mq es{p + A)+ms 

r^2 



es{p + A) + [es(p + A) + ms][eb{p) + mft] 

P , Mb - Ek(i/2) \ . 3p 



Mb + M, 



K(l/2) 



1 + 



1 E 



eb{p)+mb\ esip + A)+m, 
' IA| 



+ A) +ms 



Mb - Mk(i/2) 
>V'b(p), 



K(l/2) 



3 V Mb Ek {1/2) + Mx(i/2) 



2e„ 



d^p - I 
(27r)^ " ' ^ £^i^(i/2) + Mk(i/2) 



X , 



X 



e,(A) + 
^ 2e,(A) 
e,(A) - 
2e,(A)[e,(A) +m. 



(p-A, 

3(-E'i^(l/2) + Mk{\I2)) TTT- + 



2p 



-es(p + 



Mb - eb{p) - eq{p) + 



2(eq(p) + rUq) 



X 



Mb — Ek(i/2) 
es{A) + vfis 

2^? a1 e ( + 

Ek{1/2) + Mk{1/2) ' -Sftr(l/2) + Mif(i/2) 

^ Mb + Mj^(1/2) ^ 



Mx(l/2) 



' ^Mb^tMk^ _ \ 

^ Mb-Mk(1/2) ) 



-es(p + 



2e„ 



Ek (1/2) + Mk{i/2) 



Mb - Mk{i/2) V 
a) -e,(p + 



Mk {1/2) 

2e„ 



Ek{1/2) + -Mi^(i/2) 



•))) 



(34) 
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s2"(o) 



65 (A) + 
-es(p + 



2e„ 



2e, 



1 



i^(l/2) 



Ek(1/2) + Mk{1/2) 

~ A 



i^(i/2) s(p £;^(i/2) + -^i<:(i/2) 
)))]}^b(p), 



3 V Mb Ek{\/2) 



Z\ f d-'p - ( 



2e. 



X 



eg(A) + 
^ 2e,(A) 
eg(A) - 
2e,(A)[e,(A) + m,] 
2e„ 



— 3(-E'/f(i/2) + Mi^(i/2))- 



Ek{i/2) + ^i<:(i/2) 
p ■ A) 2p 



pA2 eq{p)+mq 



-es(p 



+ 



^+(0) = ^«(0) + (1 - ^(0) + e^f ^(0), 



Ek{1/2) + Mk(1/2) 
(2)V, 



a) -e,(p + 



es(A) + m, ^ 
2e„ 



Er {1/2) + ^ii'(l/2) 



^s(p), 



1 



^/3V Ms +M;^* 7 (2 



2e„ 



X, 



e5(p + A) 



X 



X 



\ 2e,(p + A) ^ 2e6(p) 

Mb - Ekt, 



eb{p) + rrib 



-3{Ek* + Mk, 



(P-A) 
pA2 



1 + 
1 + 



es(p + A) + ms 
Mb — Ek* 



p 



p 



eq{p)+mq es{p + A)+ms 
^2 



1 E, 



es{p + A) +ms [es{p + A) + ms][eb{p) +mb]_ 
d?p 



K. 



M, 



K* 



\ Mb Ek*+Mk: J (27r) 



* J (27r)3 V Ek* 



2eq 

+ Mk* 



(p -A) Mb- eb{p) - e,(p) 



X 



X 



eg(A) + iris 
\ 2e,(A) 
1 

2[e,(A) +m, 

^ Ek* + Mk* ' ^ Ek* + Mk* ' 



p 



3{Ek* + Mk*, pA2 2[e,(A)+m,] eq{p) + ruq 
{Mb + Mk*)[Mb - Ship) - eq{p)] + {Ek* + Mk*) 

'^b{p), 



1 Ej, 



Mk, 



Mb Ek* +Mk* J (27r)3 

-3{Ek* + Mk 



4^K*(p+^^ A) 



X 



X 



e^(A) + rus 



Ek* + Mk* 
(p-A) , p 



\ 2e,(A) 
2e,(A)[e,(A) +m. 



pA^ tq{p)+mq 



Mk* - e. 



(p+ '^-^ A) 

Er.* + Mr.* J 



Ek:; + Mk* 
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e4P + 



<(p + 



2e„ 



Ek* + Mk; 



2e„ 



Ek* + Mk^ 



Mb - Ek * 
2[e,(A)+m,]2 

2e„ 



.(p + 



Ek* + Mk* 



Mb + Mk* - eb{p) - e,{p) 



>V's(p), 



(40) 



where the superscripts "(1)" and "(2)" correspond to contributions coming from Figs. 1 
and 2, S and V mean the scalar and vector potentials in Eq. (|T8|) , ipK**,B are radial parts 
of the wave functions. Since Mk{i/2) and Mk{3/2) almost coincide with the physical axial- 
vector meson masses M„(*) and Mki we use the latter for numerical calculations. The recoil 
momentum and the energy of the K** meson are given by 



M| - Mi. 



2M, 



Ml + Ml,. 



B 



2M, 



B 



(41) 



V. RESULTS AND DISCUSSION 



We can check the consistency of our resulting formulas by taking the formal limit of h 
and s quark masses going to infinity. In this limit according to HQET the functions 



2^MbMk* 2^MbMk, VQ r , , 
^Fiw) = --j———g+{w) = — —t+{w) 



Mb + Mk* 



Mb - Mk, w + 1 



(42) 



and 



2JMbMk* 
^e{w) = — ^^s+{w), 



Mh-M, 



w 



K* 



Ml + Ml„ - e 
2MbMk*' 



(43) 



should coincide with the Isgur-Wise functions t{w) and C(^) fo^' semileptonic B decays to 
orbitally excited D mesons, B —>■ D**ev. Such semileptonic decays have been considered 
by us in Ref. |2^. Taking the formal limit — + oo, — > cxd in Eqs. (^)-(l40|) and using 
definitions (PD we find 



d^P J / 2€q 

J V^i^(3/2) (P + 



3(w + l)3/2y (27r)3 

-3M^(3/2)(u; + l)^^^ + 



M, 



K{3/2) 



w + 1) 



p 



V2 



pA^ eq{p)+mq 
d^p J / 2eg 



V'b(p) 



3 (w + 1)1/2 7 (27r)3 

-3M^(V2)(«^ + l)^^^-2 



Mk(1/2){w + 1) 



p 



eq{p)+mg_ 



i^Bip)- 



(44) 



(45) 



^As it was noted above such limit is justified only for the b quark. 
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It is easy to verify that the equahties = t and = C are implemented in our model if we 
also use the expansion in {w — l)/{w + l) {w is a scalar product of four- velocities of the initial 
and final mesons), which is small for the B — >■ D**ev decay It is important to note 
that last terms in the square brackets of the expressions for the functions iriw) (|^) and 
^e{w) ( ^51) result from the wave function transformation (|27| ) associated with the relativistic 



rotation of the light quark spin (Wigner rotation) in passing to the moving reference frame. 
These terms are numerically important and lead to the suppression of the form factor 
compared to ^p- Note that if we applied a simplified non-relativistic quark model |]7|j28[| these 



important contributions would be missing. Neglecting further the small difference between 
the wave functions 'ipK{i/2) and iIjk{2,/2)^ the following relation between and S,e would be 
obtained 



U^) = ^^^(^)- (46) 

However, we see that this relation is violated if the relativistic transformation properties of 
the wave function are taken into account. 

The relations between the form factors of heavy-to-light semileptonic and rare radiative 
B decays emerging in the large recoil limit are satisfied in our model PJT3[1. 

Using Eq. (|37D to calculate the ratio of the form factor 5'+(0) in the infinitely heavy 
h and s quark limit to the same form factor in the leading order of expansions in inverse 
powers of the heavy h quark mass and large recoil momentum |A| we find that it is equal 
to Mb I \JM'^ + ~ 0.965. The corresponding ratio of form factors of the exclusive 
rare radiative B decay to the vector K* meson -Fi(O) (see Eq. (23) of Ref. 0) is equal to 



MbI \fMB + ^K* ~ 0.986. Therefore we conclude that the form factor ratios (?+(0)/Fi(0) in 
the leading order of these expansions differ by the factor \J Mb + yM|"H^~M^ ^ 0.98. 
This is the consequence of the relativistic dynamics leading to the effective expansion in 



inverse powers of the s quark energy es(p + A) = y (p + A)^ + m^, which is high in one case 
due to the large s quark mass and in the other one due to the large recoil momentum A. 
As a result both expansions give similar final expressions in the leading order. Thus we can 
expect that the ratio r of the B branching fractions to the tensor and vector K* mesons 
in our calculations should be close to the one found in the infinitely heavy s quark limit [0 . 

The results of numerical calculations using formulas (||)-(|T0|), (|26|), (p^)-(|40|) for e = 
— 1 are given in Table |. There we also show our previous predictions for the B —* K*'y 
decay 0. Our results are confronted with other theoretical calculations and recent 

experimental data p|. The QCD sum rules predict (with 20% uncertainty) |jl9| BR{B — »• 
K*'-f) = 4.4 X 10^^ X (1 + 8%), where the second term in the brackets is the estimate of 
the terms contribution. We find a good agreement of our predictions for decay rates 

with the experiment and estimates of Ref. |^ for the measured decay rates B — > K*^ 
and B — > K2'~f. Other theoretical calculations substantially disagree with data either for 
B —>■ K*'~f or for B — > -fT^T [H decay rates. Let us note that one of the main reasons 
of the too small values for B K*^ decay rates in quark models 0,^] is the use of the 
nonrelativistic expression for the momentum of the final meson in the argument of the 
wave function overlap 0. As a result our predictions and those of Ref. [0 for the ratio 
r are well consistent with experiment, while the r estimates of and are several 
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times larger than the experimental value (see Table |). As it was argued above, it is not 
accidental that r values in our and Ref. approaches are close. The agreement of both 
predictions for branching fractions could be explained by some specific cancellation of finite 
s quark mass effects and relativistic corrections which were neglected in Ref. 0. Though 
our numerical results for the measured decay rates agree with Ref. [0, we believe that our 
analysis is more consistent and reliable. We do not use the ill-defined limit rris —>■ oo, and 
our quark model consistently takes into account main relativistic effects, for example, the 
Lorentz transformation of the wave function of the final K** meson (see Eq. (pT])). Such 
a transformation turns out to be very important and leads to the substantial reduction of 
B i^j* (1270)7 decay rate in our model. We see from Table | that our model predicts for 
the ratio BR{B K*{1270)-f)/BR{B Ki{im)-f) the value 0.7±0.3 while Ref. § gives 
for this ratio a considerably larger value ~ 2, which is the consequence of the nonrelativistic 
quark model relation (^) between form factors C,f and C,e- Thus experimental measurement 
of BR{B — > iT^ (1270)7) can discriminate between these predictions. 

VI. CONCLUSIONS 

In this paper we have investigated rare radiative B decays to orbitally excited K** mesons 
in the framework of the relativistic quark model. The large value of the recoil momentum 
|A| ~ mi)/2 makes relativistic effects to play a significant role and strongly increases the 
energy of the final meson. This effect considerably simplifies the analysis since it allows to 
make an expansion both in inverse powers of the large b quark mass and in the large recoil 
momentum of the light final meson. Such an expansion has more firm theoretical grounds 
than the previously used expansion in inverse powers of the s quark mass [0,0, which is not 
heavy enough. We carried out this expansion up to the second order and calculated resulting 
form factors in our relativistic quark model. Rare radiative B decays to axial-vector k[*^ and 
tensor K2 mesons have been considered. It was found that relativistic effects substantially 
influence decay form factors. Thus, the Wigner rotation of the light quark spin gives an 
important contribution, which leads to the suppression of the B i^i*( 1270)7 decay rate. 
In the nonrelativistic quark model, where these effects are missing, the ratio of branching 
fractions BR{B K^{1270)-^)/ BR{B fs:i(1400)7) is equal to 2, while in our model it is 
substantially smaller and equal to 0.7 ±0.3. It will be very interesting to test this conclusion 
experiment ally. 

Our predictions for the branching fractions B —>■ K*'j and B -^^27 ^^^^ their 
ratio are in a good agreement with recent CLEO data p[. 
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TABLES 



TABLE L Theoretical predictions and experimental data for the branching fractions (xlO ^) 
and their ratios Rr* = BR{B K*-f)/BR{B X,7), R ^ = BR{B k\*\)/BR{B X^^i) 

i 

{i = 1, 2), r = BR[B — > K2^)/BR[B — > K*j). Our values for the B — > K*j decay are taken from 
Ref. j. 



Value 


our 


Ref. g 


Ref. 


i 


Ref. 


i 


Ref. § 


Exp. § 


BR{B K*(892)7) 


4.5 ± 1.5 


1.35 


1.4- 


4.9 


0.5 - 


0.8 


4.71 ± 1.79 


4.55+H^ ± 0.34^* 
3.76t|]:i ± 0.28* 


Rk' (%) 


15 ±3 


4.5 


3.5 - 


12.2 


1.6- 


2.5 


16.8 ±6.4 




B (1430)7 






forbidden 










BR{B (1270)7) 


0.45 ±0.15 


1.1 


1.8 - 


4.0 


0.3 - 


1.4 


1.20 ±0.44 




Rk^ (%) 


1.5 ±0.5 


3.8 


4.5- 


10.1 


0.9- 


4.5 


4.3 ± 1.6 




BR{B i^i (1400)7) 


0.78 ±0.18 


0.7 


2.4- 


5.2 


0.1 - 


0.6 


0.58 ±0.26 




Rk, (%) 


2.6 ±0.6 


2.2 


6.0- 


13.0 


0.4- 


2.0 


2.1 ±0.9 




BR{B K| (1430)7) 


1.7 ±0.6 


1.8 


6.9- 


14.8 


0.4- 


1.0 


1.73 ±0.80 


1.66l:^J!^ ± 0.13 


Rk', (%) 


5.7 ± 1.2 


6.0 


17.3 - 


37.1 


1.3- 


3.2 


6.2 ±2.9 




r 


0.38 ± 0.08 


1.3 


3.0 - 


4.9 


0.8- 


1.3 


0.37 ±0.10 


0.39^°i^ 



* B+ K*+j 
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FIGURES 
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FIG. 1. Lowest order vertex function F*^^^ corresponding to Eq. (pi]). 
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FIG. 2. Vertex function F*^^-* corresponding to Eq. (p2|). Dashed lines represent the 
interaction operator V in Eq. (^) . Bold lines denote the negative-energy part of the 
quark propagator. 
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